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We report on the observation of weakly-bound dimers of bosonic Dysprosium with a strong universal s-wave 
halo character, associated with hroad magnetic Feshbach resonances. These states surprisingly decouple from 
the chaotic hackgound of narrow resonances, persisting across many such narrow resonances. In addition they 
show the highest reported magnetic moment ^ ~ 20 /Xb of any ultracold molecule. We analyze our findings 
using a coupled-channel theory taking into account the short range van der Waals interaction and a correction 
due to the strong dipole moment of Dysprosium. We are able to extract the scattering length as a function 
of magnetic field associated with these resonances and obtain a background scattering length = 91(16)ao- 
These results offer prospects of a tunability of the interactions in Dysprosium, which we illustrate by observing 
the saturation of three-body losses. 

PACS numbers: 34.50.-s,34.50.Cx,32.10.Dk 


Magnetic Feshbach resonances have emerged as an essen¬ 
tial tool for engineering interactions in atomic and molecular 
physics [1]. On the one hand they have been used to study 
genuine many-body systems in quantum gases, such as ultra¬ 
cold fermions in the BEC to BCS crossover [2]. On the other 
hand they have allowed major advances in few-body physics, 
through the study of the Efimov effect [3], or the production 
of ground state molecules [4, 5]. The continuous progress in 
cooling and control of atoms and molecules now opens the 
door to the study of novel condensed-matter systems using 
species with highly complex spectra. However the complexity 
of the tunable Eeshbach resonance spectrum of such species 
naturally raises the question as to whether precise control of 
interactions will be possible while avoiding harmful inelastic 
loss collisions. An illustration of this is found in the heavy, 
sub-merged shell elements Dysprosium and Erbium, both of 
which have been brought to quantum degeneracy [6-9]. Their 
complex electronic structure leads both to a large magnetic 
dipole moment and to a strong anisotropy of the van der Waals 
interaction, which give rise in turn to an extremely dense spec¬ 
trum of very narrow Eeshbach resonances, characterized by a 
chaotic distribution [10, 11]. In order to use such lanthanide 
species to study long-range many-body interacting systems, 
it is essential to understand and be able to control their colli- 
sional properties. 

In this Letter we demonstrate that broad Eeshbach reso¬ 
nances actually exist for the lowest-lying state of '®^Dy, J = 
%,mj = —8. Eurhermore these resonances decouple from the 
chaotic background and thus offer a broad tunability of inter¬ 
actions that should be readily useable for few- and many-body 
studies. We present Eeshbach spectroscopy measurements up 
to a magnetic field of 600G. This spectroscopy reveals loss 
features with widths of the order of several Gauss. Magnetic 
field modulation spectroscopy allows us to measure the bind¬ 
ing energy of weakly bound dimers. In two cases we identify 
a bound state that has a persistent universal s-wave halo char¬ 


acter over several Gauss despite crossing with a large number 
of other bound states. 

To study Eeshbach resonances we prepare ultracold sam¬ 
ples of '®''Dy atoms in their lowest Zeeman sublevel nij = —8, 
with a magnetic moment ftm = 9.93 JIb where jiQ is the Bohr 
magneton, the details of our preparation methods are pre¬ 
sented in [ 11 , 12]. The magnetic field was calibrated via radio- 
frequency (RE) spectroscopy between the two lowest Zeeman 
sub-levels of '®^Dy. Eeshbach spectroscopy is performed by 
ramping up a magnetic offset field in 15 ms to the desired tar¬ 
get value in the range from 60 G to 600 G. Typical samples 
contain 10^ atoms at a temperature of 2.4/iK (see supplemen¬ 
tal material [13] for details). The atoms are held for a wait 
time (2 s in fig. 1) in a constant optical dipole trap where they 
undergo inelastic collisions inducing losses. Eigure 1 shows 
the final atom number measured after a time-of-flight as a 
function of magnetic field with a resolution of lOOmG, map¬ 
ping the Eeshbach resonance spectrum. In this spectrum, an 
irregular pattern of several broad features appears on top of 
the sea of narrow resonances. Among them we identify in 
particular two distinct broad loss features located near 80G 
and 180G. In the following, we focus on these two features. 
A magnetic field scan with fine resolution (20 mG) over these 
two resonances shows that inside the broad features, the dense 
background of narrow resonances remain (see figure 2, and 

[13]). 

It is usually accepted that broad features in atom-loss spec¬ 
tra correspond to single broad resonances [14-17]. However, 
in the chaotic Eeshbach spectrum of lanthanides, this corre¬ 
spondence is not straightforward, and the large number of in¬ 
teraction potentials prevents so far the association of a loss 
resonance to a particular bound state, except in a few cases at 
very low magnetic field as shown by Ref. [18]. To circum¬ 
vent this shortcoming, we measure the bound states energy 
versus magnetic field using the standard technique of molec¬ 
ular association by magnetic field modulation [19, 20]. We 
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FIG. 1 . Atom-loss spectroscopy at a temperature of 2 . 4 /rK mapping the Feshbach spectrum of in 7 = 8 , ray = —8 with lOOmG 

resolution. The atom number is normalized to a reference number Aref taken at low field every 30 images. 


implement this spectroscopy on thermal clouds at a tempera¬ 
ture T with 400 nK ^ T ^ 700 nK, and focus on the regions 
B G [55G,75G], and B G [160G, 180G]. The magnetic field 
is modulated around its bias value B at radio frequencies dur¬ 
ing a time between 100 ms and 500 ms and with a modulation 
amplitude lying between lOOmG and 500 mG. 

With this spectroscopy we are able to observe the states re¬ 
sponsible for narrow resonances but we also observe an iso¬ 
lated loss resonance with typical FWHM < 30kHz. This fea¬ 
ture appears only on the low-field side of the two broad fea¬ 
tures at B ~ 80G and B ~ 180G and we associate it with 
weakly bound dimers with slowly-varying binding energy E\,. 
We extract (B) [13] and plot it in red dots in figure 2 (b). We 
put this in perspective with high-resolution atom-loss scans 
(a). In both cases we observe a very slow variation of E^ with 
magnetic field (< IMHz over lOG). This implies that these 
molecules have a magnetic moment /r very close to that of 
two free atoms, ji = 2/rat ~ ^ — 20 / 13 - This value is the 
highest reported for a magnetic dipole [18]. Furthermore, we 
observe a a quadratic behaviour for E\,{B) over several gauss, 
which reveals a coupling of the bound state with the two-atom 
continuum, conventionally termed the open channel [1]. This 
constitutes a strong evidence that the observed bound state is 
in both cases the s-wave halo state that is found in the vicin¬ 
ity of Feshbach resonances. It is striking that this strong open 
channel character is found over a broad field range, despite 
the crossing of this state with many other bound states which 
are responsible for the narrow resonances dense background. 
In fact, near the lower resonance we observe signatures of 
avoided crossings with other bound states, visible in Fig. 2, 
which correlate with narrow resonances in the atom-loss spec¬ 
trum. 

We now provide a theoretical basis supporting the existence 
of such halo states. We base the theoretical analysis (detailed 
in [13]) of these observations on the dipole-modified i-wave 
scattering model in which a strong, open-channel-dominated 
resonance overlaps with many weak, narrow resonances. Our 
observations provide clear evidence that the effect of such a 
broad resonance persists across many narrow features and that 
away from narrow poles the scattering length a takes the usual 


expression 

where Bq is the resonance pole, A its width and atg the back¬ 
ground scattering length. Furthermore in the pole vicinity, 
binding energies Bb ^ 1 MHz are well approximated by the 
universal expression: 


with m is the mass of a Dysprosium atom and h Planck’s con¬ 
stant. 

In our theoretical model we account for slight deviations 
from the universal relation (2) that result from the form of the 
long-range interaction potential which contains a typical van 
der Waals —Cg/r® term and a contribution from dipole-dipole 
interactions which for the i-wave has 0 = 1/r^ character [21]. 
We construct a two-channel model with short-range coupling 
to a ramping bound state in the closed channel. The model 
parameters are Uhg, Bq, the magnetic moment difference be¬ 
tween channel states 5/t, and the dimensionless pole strength 
■Sres abgAdjJ. [1,13]. In general, large s^es values indicate 
open-channel-dominated resonances with highly universal be¬ 
havior. 

We use this coupled-channel (CC) model to describe the 
observed data. We found out that, as expected, the observed 
resonances have a high i^es- The data allows us to extract 
only a lower bound ij-es ^10 and we cannot obtain 5/i ei¬ 
ther. Rather we can only extract the product a\,gA. We obtain 
flbg A by a direct fit to our Eh{B) data [13], the result of which 
is represented in Fig.2 (b) as a solid line. Second, we calcu¬ 
late the binding energy dependence on the scattering length 
E\,{a) assuming that abg is close to lOOao, which will be jus¬ 
tified below. We use it to translate the experimental binding 
energies to scattering lengths (red circles in Fig. 2 (c)). Com¬ 
paring the result to (1) also allows then to extract a\,gA. As 
expected these two methods agree, with a fit result given in 
Table I. In addition we also simply fit our E],{B) data by the 
universal quadratic dependence (Eb = —(B —Bo)^/(abaA)^) in 
the vicinity of the pole (B > 66G and B > 173 G for the two 
resonances respectively), with a fit result shown in dashed line 
in Fig. 2 (b) and in Table I. The coupled-channel theory and 
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FIG. 2. (color online) Broad Feshbach resonances of (a) Atom-loss spectroscopy at F = 500nK with a resolution of 20mG. This 

spectrum shows one broad feature overlapping with many narrow resonances, (b) Binding energy of weakly bound dimers measured by mag¬ 
netic field modulation spectroscopy (red circles). The solid (dashed) lines are obtained by fits of our data to the coupled-channel calculations 
(universal expression), from these fits we extract ObgA (table 1). (c) The red circles are obtained by converting the data using the 

coupled-channel calculations for a{E^). The sold lines are a fit to this data using eq. 1. The dashed lines are the scattering length resulting 
from the fit of E\,[B) by the universal expression assuming a(,g = 91 Qq (see main text). 


universal Bq 

universal ObgA 

numerical CC Bq 

numerical CC flbgA 

76.9(5) G 
178.8(6)G 

2810(100)Gflo 

2150(150)Gao 

76.8(5)G 

179.1(6)G 

2700(100) Gao 
2540(110) Gao 


TABLE I. Resonance parameters obtained from fitting the binding 
energy data to the universal formula (2) and to numerical coupled 
channels calculation. 


the universal fit are in close agreement with each other for 
the lower field resonance and in reasonable agreement for the 
higher-field one. We assume = 91 aq (as justified below) 
to calculate a{B) from the value of abgA given by the univer¬ 
sal fit, in close agreement with the CC results. Our measure 
of a is valid when the magnetic field is tuned away from the 
natTow resonances as is visible in Fig. 2 (c) (see also [13] for 
a close zoom). 

The E\,{B) data displays a quadratic field dependence well 
beyond the pole vicinity. This indicates that the observed 
broad resonances have decisive impact on the scattering prop¬ 
erties over a large magnetic field region, although multiple 
couplings to other states can be found in the data. The investi¬ 
gation of the coupling between the halo-state and other bound 
states will be the topic of future research. 

Furthermore we observe several indications that the broad 
resonance provides a background scattering length for the nar¬ 
row resonances (see [13]) that indeed crosses 0 on the high- 


field side of the 76.9G resonance. For non-magnetic atoms, 
a vanishing a results in an absence of elastic collisions at low 
temperature, rendering evaporation ineffective [22]. As a con¬ 
sequence, the final atom number and temperature after a given 
wait time are maximal. In the present case of Dy, though a 
vanishes, the dipolar interaction still induces collisions char¬ 
acterized by the dipolar length D = mpLQii^^/dinTr^ = 196ao 
and the collisional cross section = 2>2kD^/A5 [21], [? ]. 
Nonetheless when a — 0, the cross section and then thermal- 
ization is minimal. This leads to a slower evaporation, char¬ 
acterized by a maximum in the background (in-between nar¬ 
row resonances) atom number and temperature after a hold¬ 
ing time in a constant trap. We do observe a maximum 
in these two observables, located for the atom number at 
Buax,N = 109 (5) G and for temperature at Biviax.r = 107 (5) G 
(see [13] for a detailed figure). We join this observation with 
a second analysis, where we make use of the asymmetric line- 
shape of the narrow resonances. Indeed for a narrow reso¬ 
nance, the sign of its “local” background scattering length 
dbg [13, 23] influences the symmetry of its lineshape in the 
atom-loss spectrum. If dbg > 0 (dbg < 0) the local zero cross¬ 
ing attached to this resonance is on the high-(low-)field side of 
its pole, and the three-body loss spectrum is expected to show 
a minimum near such zero crossing. Most of the resonances 
we observe are too narrow to asses the symmetry of their line- 
shape given our resolution, however we can still observe sev- 
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eral ones at different fields, and they show show that below 
the 76.9G resonance, dba > 0. For higher fields we observe 
dbg < 0, up to a field of about 109 G where we actually observe 
a broad enough resonance, with a very symmetric lineshape, 
corresponding to dbg close to zero (see [13] for spectra). This 
joint analysis shows that there is an overall “background” scat¬ 
tering length that reaches 0 at a field B = 108 (5)G. Assuming 
that the broad i-wave halo resonance continues to provide the 
local B-dependent background to the narrow resonances [23], 
the universal model (1) then implies a width A =31(6) G and 
flbg =91(15) flQ for this resonance, justifying the assumptions 
made above. This value is also in agreement with the rether- 
malization measurements of ref. [24] performed at low fields. 
While the data suggests the persistence of such strong s-wave 
halo influence across a set of narrow resonances consistently 
with [23], continuing experimental and theoretical work is 
needed to explore this effect. 

We have characterized two Feshbach resonances with 
■Sres ^ 1 , we thus expect to observe effects characteristic for 
broad resonances. We now show evidence that our observa¬ 
tions of atom losses in the vicinity of the poles are indeed 
compatible with universal loss dynamics found at the center 
of broad resonances. For this we study atom-losses at differ¬ 
ent temperatures. We observe that approaching the poles, the 
final atom number reaches a minimum that is the same on the 
two broad resonances. Furthermore, at lower temperature, the 
saturation is reached in a narrower held region and at a lower 
level despite a smaller central density in the initial conditions. 
Fig. 3. 

The temperature dependence of three-body losses of Bose 
gases is known in the unitary regime where a/Ade ^ 1 , with 
AdB the thermal de Broglie wavelength. For a unitary Bose 
gas, the three-body loss parameter L 3 takes a l/T^ behaviour; 
T 3 = ^ 36-\/3 [25], [? ]. The elasticity parame¬ 

ter rjt characterizes the efficiency with which three atoms in 
contact recombine to a dimer and a free atom. 

We compare our findings with the loss dynamics model de¬ 
veloped in [26], which predicts the final atom number taking 
into account two-body evaporation and three-body recombi¬ 
nation. This model requires the knowledge of the trap depth, 
frequencies and initial atom number and temperature. These 
observables are calibrated (see [13]), and the only unknown 
parameter of this model is Tj*. We then compare this model 
with a fixed 77 * to the three data sets represented in figure 3. 
Within experimental uncertainties, we And that this model re¬ 
produces well our Anal atom number at resonance for a fixed 
rjt that we And to be 

t7*=0.07+[|;i^. (3) 

This value is close to the lowest reported values for alkali- 
metal atoms [27-32] and implies lifetimes comparable to 
these species in the unitary regime. Having this experimental 
information is valuable since no theoretical prediction is 
available. We observe that Tj* is identical for our two broad 
resonances. A systematic study of L^(T) at resonance would 



FIG. 3. (color online) Atom loss spectroscopy for different exper¬ 
imental starting conditions and wait times. Final atom number N 
normalized to the initial one Nq as a function of magnetic fleld, with 
a resolution of lOOmG. Green: initial temperature: Tq = 2.4pK, ini¬ 
tial density: hq = 6 .1(2.0) x 10*^ cm^^ wait time t = 2s, blue: Tq = 
0.5pK, no = 3.7(1.3) x lO'^cm^^^ f = 500ms, red: Tq = 0.3pK, 
no = 4.6(1.6) X 10*^cm^^ t = 300ms. The temperature dependence 
of the saturation is well reproduced by the model of universal loss 
dynamics of unitary Bose gases of ref. [26] (solid horizontal lines). 
The shaded regions represent the uncertainty on the results of the 
model given a one-standard deviation on all experimental parame¬ 
ters entering the model (temperature, atom number, trap frequency 
and depth). 

yield a more precise measure. 

In conclusion, it is remarkable that despite an apparent 
chaotic distribution of Feshbach resonances, isolated states 
can decouple from this background and lead to a simple 
description. This emergence of structure from a chaotic 
background is a generic feature of quantum-chaotic systems 
[33], with the emergent states associated with classical 
quasiperiodic trajectories in phase space [34]. Further work 
is needed to understand the origin of our broad resonances. 
However the simple universal description provided here 
is sufflcient to predict how to control dipolar gases in an 
understood way [35]. Between narrow resonances we obtain 
pure BECs with up to 25 x 10^ atoms with long lifetimes 
(> 1 s). The scattering length obtained in this work can 
thus be conflrmed with precision by the measurement of the 
density distribution of Bose-Einstein condensates. Such mea¬ 
surements are ongoing with our setup. This opens prospects 
of studying the few-body physics of bosonic dipoles; for 
instance the energy of Eflmov states is thought to be impacted 
by the dipole-dipole interaction [36]. 
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Supplemental Material 

EXPERIMENTS 
Atom-loss spectroscopy 

The early stages of our experiment have been presented 
in [ 1 , 2 ], we detail here the specifics of our atom-loss spec¬ 
troscopy. The atom-loss spectroscopy consists in recording 
the number of atom left after a given wait time as a func¬ 
tion of magnetic field. On a Feshbach resonance, three body 
losses are amplified so that it shows as a minimum in final 
atom number. The parameters used for this technique vary 
in the different data sets presented throughout the main text. 
All experiments are performed in a crossed optical dipole trap 
created by a narrow-line laser with wavelength 1064nm. We 
obtain the trapping potential created by this laser using stan¬ 
dard techniques. A forced evaporation is performed in this 
trap down to the temperature of choice by lowering the beams 
power. The final power is varied to control the temperature T 
of the gas. We express the trap depth U through the trap ratio 
r] =U/k^T. It can be estimated using the knowledge of the 
beam waists and the effect of gravity. This is strongly depen¬ 
dent on the knowledge of the trapping potential, furthermore 
since dysprosium is strongly magnetic, residual magnetic field 
gradients can change the depth. We then improved our estima¬ 
tions using a statistical analysis of final temperature and atom 
number in our atom-loss spectroscopy using loss dynamics 
equations adapted from [3] for finite s-wave scattering length 
and including universal dipolar scattering. We obtain a knowl¬ 
edge of 77 with an uncertainty of 20 %. 

In figure [1] of main text we present a spectroscopy per¬ 
formed at a temperature of 2.4/rK. The samples were pre¬ 
pared in the following way; after evaporation at a magnetic 
field close to OG down to T ~ 0.6 fiK, we recompress the trap 
by increasing the beams power by a factor of approximately 
x2 in 100 ms. This has two purposes, first we reduce resid¬ 
ual evaporation induced by increased elastic collisions, second 
we minimize a weakening of the trap confinement caused by a 
residual field gradient created by our Feshbach coils, and par¬ 
tially to the crossing of many Feshbach resonances. We then 
increase the bias magnetic field in 15 ms to the target value. 
We observe that the turning-on of the magnetic held to val¬ 
ues above 60 G results in a residual heating, this heating is 
partly due to an unintentional change of direction of the mag¬ 
netic held when ramping to high helds. At the change of sign 
the cloud is in a zero magnetic held and is slightly depolar¬ 
ized which results in dipolar relaxation and heating when the 
held is increased. A residual magnetic held gradient caused 
by our Feshbach coils probably contributes to this heating as 
well. The initial temperature of the clouds for the atom-loss 
spectroscopy is T = 2.4/rK, and the starting atom number 
No = 1.8 X 10^. We estimate the trap depth in these condi¬ 


tions to be rjo = k^To/U = 7.0(1.5). The hnal atom number 
presented in hgure [ 1 ] of main text is recorded after a wait 
time of 2 s at constant magnetic held. 

In hgure [2] of main text, the atom-loss spectra presented 
were obtained at a temperature of 500 nK with initial atom 
number Aq — 1.0 x 10 ^, the two data sets are taken at a similar 
trap depth of t] ~ 5. The wait time is in both cases 500ms. The 
data set spanning the lower resonance was taken with a higher 
mean trapping frequency than the other one which induces a 
higher initial density, thus this data shows a lower hnal atom 
number. 

In hgure [3] of main text, we present three different data 
sets, in different experimental conditions. The data presented 
in green (color online) is the same as in hgure [ 1 ] of main 
text. The parameters for the three different data sets presented 
(green, blue and red) are summarized in table II. We use these 
parameters to solve numerically the coupled equations ( 12 ) 
and (25) of [3] from which one obtains the expected hnal atom 
number for unitarity-limited two and three-body losses in a 
Bose gas. These equations read 


dT 

dt 


T 

3 


( 5)6 
V 3 rs 


dN „ Vev 


(4) 

(5) 


with 


nk^U 


72 =h 


mO) 


2\/37tk^ 


A3 = Li 


(6) 

(7) 

( 8 ) 


where L 3 is given in the main text and the expressions of ^ 
and k can be found in ref. [3]. These equations can then easily 
be solved numerically remembering that 77 = T/U, using the 
parameters given in table II. 

These predictions are compared to our measurements in hg¬ 
ure [3]. The shaded area represents the range of calculated 
hnal atom number we obtain by varying the input parameters 
of the model by one standard deviation as given in tab. II. 


Localization of the zero-crossing 

We here employ the (a)symmetry of narrow Feshbach res¬ 
onances lineshapes to extract their local background scatter¬ 
ing length, in the surroundings of the 76.9G resonance, rep¬ 
resented in hgure 4. The sign of the local background scatter¬ 
ing length on a narrow resonance dbg determines the side on 
which the local zero-crossing takes place, if dbg > 0 (dbg < 0 ) 
it is on the high- (low-) held side of the resonance. The zero¬ 
crossing corresponds to an attached minimum in three-body 
recombination in its vicinity, that appears as a maximum of 
atom number in the atom-loss spectra. If the resonance is suf- 
hciently broad to extend over several data points, this leads to 
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FIG. 4. Atom-loss spectrum in the surrondings of the broad resonance at 76.9 G. We observe an inversion in lineshape when crossing this 
resonance (between the two first panels), and then a symmetric lineshape at 109 G suggesting that the local background scattering length 
changes sign in this region. 
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FIG. 5. Zoom in on narrow resonances, (a) and (b) show resonances with a zero-crossing on their high-field side (5t,g > 0) which occurs for 
fields B < Bq = 76.9 G. (c) and (d) are examples with inverted asymmetry occurring for B > Bq, implying dt,g < (note that the resonance in 
(e) overlaps with several very narrow lines). Finally at 109G (e) we observe a symmetric lineshape which corresponds to a a^g close to zero. 


an asymmetry of its lineshape in the data, as presented in fig¬ 
ures 4, 5 also visible for instance in spectra from [4]. We are 
able to extract the sign of dbg on several resonances at differ¬ 
ent magnetic fields, and observe the following: For fields up 
to the position of the broad resonance at 76.9G, all the line- 
shapes of isolated resonances we can observe correspond to 
a dbg > 0 (fig. 5, (a, b)). The symmetry of the lineshapes is 
inverted when crossing the broad resonance, as seen in fig.5 
(c,d). Finally, at a field B = 109 G we observe a resonance 
with a very symmetric lineshape, which suggest a dbg close to 
zero, we conclude that a change of sign of dbg occurs in this 
field region. We do not observe sufficiently broad spectra be¬ 
tween HOG and HOG. Higher fields are closer to the other 
broad resonance at 178.7G. Since it is a broad resonance, 
it likely provides the background for the narrow ones in its 
vicinity, and we cannot conservatively infer information about 
the 76.9G resonance in this region. Our observation of a zero 
crossing of the background scattering length near B — 109 G 
is corroborated by background atom number and tempera¬ 


ture measurements. We present in figure 6 the atom number 
and temperature as a function of magnetic field in an atom- 
loss measurement. We observe that the background value 
of these observables (in between narrow resonances) varies 
slowly with magnetic field and reaches a maximum. From 
the data in fig. 6, we locate the maximum in atom number 
at BMax.iV = 109 (5)G and temperature at Bjviax.r = 107 (5)G. 
This maximum in atom number and temperature signals a 
minimum of thermalization due to a minimum of collision 
rate, which hampers evaporation. The two-body collision rate 
is related to the elastic scattering cross-section, which is for 
dysprosium the sum of the constant universal dipolar scatter¬ 
ing cross section and that of the short-range s-wave collisions: 
fftot = 327rD^/45 + where D is the dipole length given 
in main text. Since the dipolar scattering cross section is in¬ 
dependent of magnetic field, the observed minimum in col¬ 
lision rate corresponds to a minimum of a^. This minimum 
appears to be located close to the field at which the local back¬ 
ground scattering length is close to zero according to the line- 
















8 



95 100 105 110 115 


B(G) 

FIG. 6. Focus on the background atom number and temperature in 
between narrow resonances after a wait of 500 ms, A^ref and Tj-ef are 
taken after the same wait time at a field close to zero. A maximum 
is visible on both observables, that we locate for the atom number 
at BMax.N = 109 (5) G and for temperature at fiMax.r = 107 (5)G, 
shown as grey areas. 


data set 

Ao(10") 

7b(nK) 

Vx,y,z (Hz) 

m 

fwait (ms) 

green 

1 .8(2) 

2400 (200) 

70 (10) 
123(15) 
263(25) 

7.0(1.5) 

2000 

red 

0.46(2) 

500(100) 

44(5) 

87 (10) 
145(15) 

4.0 (1.0) 

500 

blue 

1 .0(2) 

450(100) 

52(5) 

23(5) 

260(25) 

5.0(1.0) 

500 


TABLE II. Experimental conditions for the atom-loss specta pre¬ 
sented in figure [3] of main text. Nq and Tq are the initial atom 
number and temperature, rjo = U/k^To where U is the trap depth. 
Vxj^z refers to the trapping frequencies in the three dimensions, fi¬ 
nally t^ait is the time we wait before recording the final atom number 
N 

shapes analysis. With this combined analysis, we conclude 
that indeed the scattering length in between narrow resonances 
varies slowly and crosses zero at a field B = 108 (5) G. 


We associate very broad spectra to bound states with high 
energy-vs.-magnetic field slope, meaning that their magnetic 
moment differs strongly from that of two free atoms. On the 
opposite, the center frequency of the narrow feature evolves 
slowly with magnetic field and we associate it with the weakly 
bound state attached to our broad resonances. We present in 
figure 7 (right) two spectra on each resonance from which 
we extract the binding energy using a fit with a Maxwell- 
Boltzmann distribution, solid lines in figure 7. This fit ac¬ 
counts for the broadening of the spectrum by the relative ki¬ 
netic of two atoms associated. The fit results are typically of 
good quality given our noise level, the width observed are in 
qualitative agreement with our temperature. A residual heat¬ 
ing of about 20% increase in T is observed when dimers are 
associated, which we attribute to energy released by the relax¬ 
ation of shallow dimers to deeply bound states. This dynam¬ 
ics surely modifies our spectra, but since the binding energy 
is given by the activation energy, meaning the lower bound of 
the line, its measurement is independent of details of the line- 
shape. To account for possible systematic errors, we still put 
conservative error bars on our data. Finally, in this figure one 
can observe a weak association feature on the left panel at low 
frequencies, this is the sub-harmonic of the main association 
line. 

Avoided crossing regions 

In the binding energy data we observe avoided crossings 
that modify locally the halo-dimer binding energy. They result 
from a stronger-than-usual coupling between the halo dimer 
and another state. These influence also the scattering length 
and we thus also observe their effect in the atom-loss spec¬ 
trum. In these regions the scattering length obtained in the 
main text does not hold true, outside of them all resonances 
are very narrow (< 20 mG in the atom-loss spectrum) and per¬ 
turb the scattering length only very locally. In figure 8 we 
outline in blue the regions where the extracted a{B) does not 
hold. 


Binding energy measurements 

We measure the binding energy of weakly bound dimers 
via a modulation of the magnetic field around its bias value. 
This modulation is applied with an additional coil placed near 
the atom cloud, with a current oscillating at radio-frequencies. 
We record the number of atoms left after the field modulation 
as a function of frequency at fixed bias field. The typical am¬ 
plitude of modulation is 300 mG during a time of 300 ms. The 
output of this method presents typically two kinds of spectra, 
which are exemplified in figure 7 (left). We observe either 
losses in a narrow frequency band, or in a very broad band. 
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FIG. 7. Magnetic-field Modulation spectroscopy. Left: output of magnetic-field modulation, in the magnetic field range B G [72G, 72.6G]. 
Number of atoms normalized to the number without field modulation Nq, as a function of magnetic field and modulation frequency. We observe 
one narrow, slowly varying feature and one localized very broad feature. The former is the signature of the weakly bound state creating the 
broad Feshbach resonance at 76.9 G. The latter is created by a bound state responsible for a narrow Feshbach resonance. Right: Examples 
of magnetic-field modulation spectra from which we obtain the binding energy £ 1 , shown in figure [2] of main text. The data is fitted by a 
Maxwell-Boltzmann distribution to extract E^,. 



FIG. 8. Zoom on the 77G resonance, same as fig. [2] of main text. 
We outline here in blue the regions where we observe strong avoided 
crossings in the binding energy as well as stronger-than-usual nar¬ 
row Feshbach resonances. In these regions the scattering length is 
unknown. Outside of them resonances are very narrow and modify a 
only very locally. 


THEORY 

Long-range interactions of Dy atoms 

Collisions of ultracold lanthanide atoms such as Dy are 
very challenging to describe theoretically due to the open 
4/ electronic shell and the resulting large electronic angu¬ 
lar orbital momentum. Another difference between Dy and 
alkali metal atoms is the strong magnetic dipole of 9.93/fB- 
The anisotropic dipolar and van der Waals interactions are so 


strong that “partial wave” L is not a good quantum number 
for the eigenstates of Dy 2 molecules, but many different elec¬ 
tronic symmetries and partial waves are mixed by the strong 
anisotropy. Furthermore, the density of states is very high, 
corresponding to a mean spacing between narrow resonances 
on the order of 0.3G for the bosonic species of Dy [2, 5] as 
well as the similar species Er [4]. In fact, the series of very 
narrow resonances observed by Frisch et al. [4] is known to 
have a mean spacing explained by the Wigner-Dyson distribu¬ 
tion predicted by random matrix theory, implying a “chaotic” 
set of levels that do not have any assignable set of quantum 
numbers. Therefore, the appearance of the broad identifiable 
i-wave bound states in the present experiment is somewhat 
of a surprise, since this well-characterized universal state is 
quite different in character from the set of “chaotic” levels 
that make the dense set of narrow resonances. The long-range 
s-wave “halo” state effectively decouples from the “sea” of 
chaotic levels and generates a broad resonance that spans a 
relatively large tuning range of B. In this case, for the mag¬ 
netic field regions of our current interest, the impact of narrow 
resonances is negligible outside very small regions around 
their locations. Each wide resonance provides a “local” back¬ 
ground for nearby narrow ones, as illustrated for alkali atoms 
in Refs. [6,7]. 

The experiment prepares the Dy atoms in the maximal pro¬ 
jection state with angular momentum 7 = 8 and projection 
m = — 8 . The dipolar physics of two Dy atoms can be set 
up in different basis sets. Kotochigova and Petrov [ 8 ] set 
up the usual body-frame Born-Oppenheimer electronic states, 
whereas Bohn et al. [9] set up a space-frame asymptotic ba¬ 
sis in a partial wave expansion. Since we are primarily con- 
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cerned here with the long-range properties of universal s-wave 
states, we will use the latter to describe the collisions of two 
Dy atoms in the same Zeeman sublevel. This requires a single 
i-wave entrance channel, although we will need to consider its 
off-diagonal coupling to the t/-wave channels. The many other 
electronic and partial wave channels only come into play by 
generating a dense set of threshold crossing levels with very 
different magnetic moments that are weakly coupled to the 
s-wave halo state. 

The basic understanding of the s-wave halo state requires 
only a knowledge of the long-range potential, which is com¬ 
prised of one or more terms of the form Cp/R^. Here p can be 
6 , 3, or 4 for the respective van der Waals, dipolar, and non- 
adiabatic contributions [9] . There are several definitions in the 
literature for the scale length of a CpjRP long-range potential, 
and we adopt here the definition of Gao [10] 


Rp = 




(9) 


Note that the lengths /?vdw d defined for the p = 6 van 
der Waals potential in Ref [11] are respectively jRe and 
0.477989R6. Note also that the dipole length D defined by 
Bohn et al. [9] is ^R^. These lengths define corresponding 
energies Ep as in Eq. 9. 

Bohn et al. [9] consider the collisions of two dipoles with 
their orientation held fixed in space by an external field. In our 
case, each atom maintains its space projection mj = —8 when 
the distance between two colliding atoms is large enough that 
the Zeeman splitting due to the laboratory B field is large com¬ 
pared to the diagonal and off-diagonal potential terms due to 
the van der Waals and dipolar interactions. This is the case 
at distances larger than the shortest length scale Rg needed 
or analyzing the i-wave halo state. Thus, each atom stays in 
the same \j,m) = | 8 , — 8 )projection state and the asymptotic 
basis set is the set of all states | 8 , — 8 )i | 8 , — 8 ) 2 |TMi), which 
we need to label only by the partial wave quantum numbers 
L,Mi. For identical bosons we need only even L values, and 
only partial waves with Ml = 0 are in the same Hamiltonian 
that contains the s-wave, since the total projection quantum 
number —16 is conserved. Treating the short-range states and 
the dense set of chaotic eigenstates would require including 
the full basis set, \jmi)i\jm 2 ) 2 \LML) for all separated atom 
combinations and the mixing of many partial waves [4], but 
we can ignore this complexity in estimating the effect of s- 
wave universality, and incorporate the effect of such complex 
states as causing the narrow avoided crossings with the uni¬ 
versal i-wave halo state in the near-threshold region. 

The asymptotic space-frame physics is described by Eqs. 
(2)-(5) of Bohn et al. [9], where we will estimate the ef¬ 
fect of dipolar interactions using only the \LMl) = |0,0) and 
12,0) terms in the Hamiltonian, described by the coupling 
coefficients (in the notation of Ref. [9]) = 0, Cq 2 ^ = 

l/>/5, and = Ijl. Thus, the diagonal dipole interaction 
—C 3 /R^vanishes for the s-wave, and is C 3 = foJ" 

diagonal L,Ml =2,0 (7-wave component (the other (7-wave 


components have ^ 02 ^'^ ~ —Ijl and con¬ 

tribute to universal dipolar scattering, but do not couple to 
the s-wave). Since for Dy d = 9.93/J.b and ps = Ct/2 atomic 
units, where a is the fine structure constant, we find that 
(7^ = 0.00131 atomic units (1 atomic unit = where Eg 
is one Hartree and ao is the Bohr radius). The Bohn scaled 
unit of length D — 196ao and the Gao one is r-j, = 392 ag for 
C 3 = (71 

Although dipolar interactions vanish in the diagonal s-wave 
potential, Bohn et al [9] show how a simple perturbation es¬ 
timate shows that the off-diagonal coupling between the s- 
and (7-wave Mi = 0 components give an effective adiabatic 
potential that varies at long-range as —Ci,lR^. This is be¬ 
cause the splitting between these components at long-range is 
given by the centrifugal potential, 6P/{2p.R^), which is much 
larger than the off-diagonal dipolar term {2/p5)d'^/R^. The 
long-range adiabatic coefficient is C 4 = (4/15)(7^Z) = 0.0687 
atomic units E^a^. This gives a Gao scale length of R 4 = 
143 ao. This is comparable to the scale length of the van der 
Waals potential, for which Ce = 1890 au [ 8 ], Rg = 154 ao, and 
E(,/h = 0.9265 MHz. 


Overlapping Feshbach resonances 

Many alkali metal atom systems are known to have broad 
i-wave resonances on which are superimposed one or more 
much narrower resonances of higher partial wave bound state 
character; for example. Ref. [ 6 ] for Cs -f Cs gives an example 
with an overall very large aig and Ref. [12] for Rb -f Cs gives 
an example with a quite small a^g. Reference [7] developed 
an analytic formalism based on multichannel quantum defect 
theory (mqdt) to characterize such sets of overlapping reso¬ 
nances. The mqdt treatment assumes the following: (1) one 
i-wave entrance channel with “background” scattering length 
ahg and p — 6 length Rg, (2) a set of A > 1 closed channels 
each with an uncoupled “bare state” with crossing positions R, 
ramping according to magnetic moment differences 5/r, from 
the entrance channel, and (3) coupling between the entrance 
channel and closed channel i described by a dimensionless 
coupling matrix elements y,, or equivalently, a set of “bare” 
resonance “widths” A,. The result of the analysis is that the 
scattering length with a set of N poles can be put in the fol¬ 
lowing form: 

a{B)^aigfl(^l-j^y (10) 

Here atg provides a “global” background for all the reso¬ 
nances, and the “widths” A, and resonance pole positions RT® 
are generally functions of all the “bare” widths A, and crossing 
positions R,. 

One consequence of the general formula is that near any 
particular pole i we can write 

/X / A; 

a(R) = dbg,i [1 - g gres 


( 11 ) 
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where the “local” background scattering length is 

abg,i = “fog 0 ~ gres ~_gres ^ 

Here only the product UbgjAi is well-defined, whereas the in¬ 
dividual terms can depend on the number of resonances and 
B field range included in the model. Thus, the “local” pole 
strength of a resonance is well-defined, and certainly there are 
cases where a broad i-wave resonance can be characterized 
by a single ig parameter as B is tuned across several narrow 
resonances. A good example is provided by the pair of s- 
and g-wave resonances of two Cs atoms with respective po¬ 
sitions and pole strengths of {Bi,Sres.i) = (548.44G, 160) and 
(554.060,1.6) [7]. Accounting for the overlapping proper¬ 
ties of these two resonances is essential in characterizing the 
three-body physics in the 550 G region of Cs [13]. 

Figure 9 gives an example of how the scattering length can 
look within this treatment for the case of a single very broad 
resonance overlapping with many narrow ones. In this exam¬ 
ple with a toy model, a random distribution of narrow reso¬ 
nances was taken with two different values assumed for the 
mean width. The broad resonance provides the “local” back¬ 
ground dbgj for the narrow resonances. The range over which 
a narrow resonance permrbs the underlying background due to 
the strong resonance depends on the width of the narrow one. 
Figure 8 indicates regions where such perturbations are evi¬ 
dent in the binding energy curve. It will be important for fu¬ 
ture theoretical and experimental work to investigate the pole 
strengths of such narrow resonances and assess their effect. 
Reference [14] describes how it is not valid to extrapolate the 
properties of a broad resonance across a range on the order 
of the spacing between different vibrational levels of the same 
channel, although it is valid to extrapolate its properties across 
a set of resonances in different spin channels. 

Universal long-range theory for hybrid interaction potential 

The scattering length a{B) and the energy E{B) of the last 
i-wave bound state near an isolated magnetically tunable near¬ 
threshold Feshbach resonance is given by [11]: 

0 = -^ 

where Obg is the background scattering length of the s-wave 
entrance channel in the absence of the resonance, B is mag¬ 
netic field strength, Bq is the pole position of the resonance 
as collision energy E ^0, and A represents the width of the 
resonance. The “bare” closed channel level that makes the 
resonance ramps linearly with B as 5li{B — Be), where 5/i is 
the magnetic moment difference between the closed channel 
state and the two separated atoms. The actual pole position 
Bq is shifted from the “bare” crossing location Be by a shift 
that can be as large as A, depending on Obg [11]. While the 



FIG. 9. Scattering length for the case of a single resonance similar to 
the 77G one observed in the experiment overlapping with 50 much 
narrower ones (upper panel: mean width of narrow resonances is 
0.125 G, lower panel: mean A equals 0.004 G). The red thick curve 
shows the background induced by the wide resonance. 

formula for a{B) applies over a wide range of B, the “univer¬ 
sal” binding energy formula for E{B) only applies very close 
to threshold where a ^ R(,, where Rg is the scale length of the 
van der Waals potential. 

The “pole strength” of resonance term in Eq. 13 is charac¬ 
terized by the product abgA and more generally by a dimen¬ 
sionless parameter sg = (agg/Rg)(5/rA/£’g), or equivalently, 
by the s^es = 0.477989sg parameter used in Ref. [11] to clas¬ 
sify resonances (see also Ref. [14]). Strong or “broad” reso¬ 
nances have ires ^ 1 and weak or “narrow” resonances have 
■Sres ^ 1- The former have universal i-wave bound states as B 
is tuned over a large fraction of the resonance width A whereas 
the latter are universal only over a small fraction of their (gen¬ 
erally small) A near the pole. Reference [15] studied in detail 
the departure from universality of the binding energy for broad 
and narrow resonances in the ®Li and ^Li systems. 

Using the Gao scaled units, r = R/Rp and £ = EjEp, the 
single channel Schrodinger equation for a power law —Cp/R^ 
potential becomes 

^^-Vsr{r)^W = Q, (14) 

where Vsr{r represents the short-range contribution. For nu- 


^ / 1 _L(L 

dr^ \ rP r 
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FIG. 10. (Upper Panel) Single channel binding energy curves E jEp 
versus scattering length a/Rp for p = 4, 6, and for a hybrid potential 
with both p = 4 and 6 terms to simulate the i-wave collision of two 
weakly dipolar Dy atoms. The C 4 /C 6 ratio for the latter is chosen 
to be the value from the perturbation analysis for two Dy atoms (see 
text). The hybrid case is scaled by Ug and R(, from the p = 6 term. 
The dashed line shows the bottom of the “bin” for the last bound state 
of the van der Waals potential at E/E(, = —157.8. (Lower Panel) 
Same as above but on a log scale near threshold. The dashed line 
shows the universal quadratic relation in Eq. 13, to which the differ¬ 
ent curves all tend near threshold when a/Rp > 10. The departures 
from quadratic universality follow similar curves for the p = 4, 6, 
and hybrid potentials for a/Rp > 3, but these cases show variation 
from one another at smaller a. The crosses show the coupled chan¬ 
nels curve from the 2-channel calculation with Jres = 16. For Dy, 
Ee/h = 0.926 MHz and R^ = 154 ao. 


merical purposes we can simply represent this short-range part 
as a Lennard-Jones p — 2p potential, as Vsr{R) = +(7^1 
where the short-range parameter a is adjusted to select the 
number of bound states in the potential and the scattering 
length. When the number of bound states varies from n to 
n-F 1 as (7 is varied over a range (7„ to OF+i, the scattering 
length of the power law potential can be varied from —00 to 



FIG. 11. Plot of the information in the lower panel of Fig. 10 in the 
form Ea^/{E(iR^), which approaches unity for the universal relation 
in Eq. 13. The clear departure from the universal scaling due to the 
finite range of the potential is evident in the range of a/Rf, < 10, 
which is the range sampled in the Dy experiment, indicated by the 
data points, where each point at 6 = Bi is plotted using the a(B,) 
mapping given by the coupled channels calculation. The departure of 
the data from the predictions of the hybrid potential model at smaller 
fl/Bg is likely due to the breakdown of our over-simplified 2-channel 
model as binding energy becomes larger and a smaller. 


4 - 00 . Similarly, the binding energy of the last bound state can 
be varied from its largest possible value at the bottom of the 
last “bin” in which an s-wave bound state must be found to 
zero energy when the scattering length approaches 4-°°. Such 
a calculation gives the E{a) dependence for a single chan¬ 
nel and describes open-channel-dominated resonances with 
■Sres 7^ 1. In this sense, open-channel-domianted resonances 
are all universal, although they show potential-dependent de¬ 
viations from the qudratic universality described by (13). 

Figure 10 shows the Sres ^ 1 binding energy curves 
E{a/Rp)/Ep for p = 6 and p = 4 in scaled units. These curves 
are “universal” in that they apply to any species and are in¬ 
dependent of reduced mass and the C„ coefficient by using 
scaled units. The bottom of the “bin” for the last bound state 
is at — ISSfig and —105£’4 for the respective p = 6 and 4 po¬ 
tentials. The Figure also shows the E{a) curve, scaled by Eg 
and ag, for the hybrid adiabatic i-wave potential for Dy with 
Cg = 1890 au and C 4 = 0.0686 au, as described previously. In 
this case, E(,/h = 0.926 MHz and /?g = 154 ap. 

The lower panel of Fig. 10 shows an expanded view on 
a log scale of the region near threshold, where the conven¬ 
tional “universal” quadratic binding energy curve applies in 
the limit a ^ based on the approximation in Eq. 13 that 
E/Ep = (Rp/a)^. Figure 11 shows the same information plot¬ 
ted as Ea^/{E(iR^) versus a/R^, which approached unity as 
a —> 00 (see Ref. [15]). The actual binding energy curves are 
very close to the universal form when the scattering length is 
IQRp or larger, but begin to show noticeable variation from the 


































13 


universal curve when the scattering length is smaller. The ac¬ 
tual binding energy differs by a factor of 2 from the universal 
curve when the scattering length is as small as 2R(,. Therefore, 
a much better mapping of measured binding energy to scatter¬ 
ing length will be obtained if the actual scaled E{a) curve is 
used instead of the universal one. 


Coupled channels model of broad Feshbach resonances 


A full coupled channel model describing the properties of 
the wide resonant features observed in the experiment is be¬ 
yond our reach. Instead, we turn to a simplified two-channel 
description, neglecting the impact of the natTow resonances, 
which in this case seems to be a good assumption based on the 
experimental data, namely, that a universal i-wave halo state 
emerges from the sea of chaotic levels due to the ramping with 
B of a particular approximate eigenstate to which the entrance 
channel is strongly coupled. We construct a 2-channel model 
using the method described in Refs. [13, 16, 17] for reducing 
the full coupled channels model with many resonances to an 
effective 2-channel one based on potentials having the right 
long-range form, but with a reduced number of bound states. 
We assume a Hamiltonian matrix of the form; 


( V,{R) Arese-^’^ \ 
\Arese-^^ Vc{R)+Eres ) ’ 


(15) 


where Vg{R) describes the hybrid ground state adiabatic po¬ 
tential with the 6-12 LJ form plus the p — A term, the off- 
diagonal term describes the coupling between the entrance 
and the closed channel with potential Vc{R) -f and E^es > 0 
is the asymptotic energy of the closed channel; we assume Vg 
and Vc both vanish as B oo. The inner wall of Vg is varied 
to select the scattering length and the number of bound states 
in the entrance channel. We chose to set ai,g = 102 ao and 
the number of bound states to be 13. As in Refs. [16, 17], 
we take Vc to be the same potential as Vg. The —2 vibrational 
level with a binding energy of 230.507 MHz was selected to 
represent the Feshbach level; thus, the “bare” Feshbach level 
will be at the E — 0 threshold of the entrance channel when 
Eres/h = EccJh = 230.507 MHz. Varying Eres thus allows us 
to “ramp” the bare Feshbach level across threshold and tune 
the scattering length of the coupled system, which will gener¬ 
ally have a pole at a shifted position E^cs ^ assume 

Eyes is a linear function of B with slope 5/^, we can select the 
mapping of experimental B to the Hamiltonian model param¬ 
eterized by Eyes'. 

Eyes = E‘yes + 5p{B- Be) = E^s + 5p{B- Bo) (16) 


where By is the field where the bare Feshbach level crosses 
threshold. The pole of the scattering length generally occurs 
at a different field Bq, with the shift from By being described 
in Ref. [11]. 

The inter-channel coupling V],c(^)is characterized by the 
off-diagonal term proportional to Ayes, which is selected to 


give the resonance “pole strength,” or Sres parameter; here we 
take j3 = 1 ag *. The energy-dependent width HY{E) of the res¬ 
onance for collision energy £ > 0 is (see Eq. 36 of Ref. [11]) 

hr{E) =27t\{b\Vgy{R)\E)\^ = 2{ka){Esyes) (17) 


where \b) and \E) are the wave functions of the respective 
bound state and energy-normalized s-wave scattering state, hk 
is the relative collision momentum; equivalently, we could de¬ 
fine fiY{E) = 2ka(,EoSf,. 

Given that we have fixed a value of a^g, which we here 
take to be “nominal” at a magnitude order of 100 ao, the free 
parameters in the model are Eyes for resonant tuning and the 
strength of off-diagonal coupling Ayes, which selects the “pole 
strength” ires- Equation 16 can be used to map the tuning 
parameter to a laboratory B field, given a value of 5p and 
Bq. The ires parameter determined by selecting Ayes is also 
proportional to 5p, since by definition (Eq. 35 of Ref. [11]) 


abgASp 

aE 


(18) 


The 2-channel coupled Schrodinger equation is solved numer¬ 
ically using standard methods and the bound state energies for 
E <0 and i-wave scattering length for £ 0 are easily found 

as a function of Eyes and ires (or equivalently. Ayes), given the 
fixed djyg. 

Eigure 12 shows examples of £ versus Eyes for different ires 
parameters. The upper panel shows the bare states (uncou¬ 
pled, Ayes = 0) and the coupled eigenstates for ires = 8.719. 
The avoided crossing and the shift in pole from the bare cross¬ 
ing at Eyes/h =230.51 MHz to the pole at Eyes/h = 224.14 
MHz are evident. The lower panel shows the binding energy 
of the threshold bound state versus Eyes — Eyes, that is, relative 
to its pole position, for several different values of ires- The 
width of the resonance increases linearly with ires (by defini¬ 
tion from Eq. 18). In the ires —> 0 limit, the binding energy 
curve would approach the linear curve for the bare state, with 
a small quadratic departure only very near threshold. Keep in 
mind that the near-threshold £ versus a curve becomes univer¬ 
sal, that is, independent of ires when ires ^ T and approaches 
the single channel curve, as indicated in Eig. 10. The univer¬ 
sal £ (a) curve is independent of 5p when ires ^ 1, in practice 
where a/Ro > 1 and ires ^10. Consequently, a is determined 
by measuring £ in this universal limit (away from perturba¬ 
tions by narrow levels, of course). 

Einally, Eig. 13 shows a comparison of a coupled chan¬ 
nels calculation with the measured binding energies versus B 
for the two resonance regions. The fit is not unique, since 
ires and dp are not determined independently for such a halo 
state. The figure shows fits with different assumed 5p and 
ires values. Taking 5p — gps with g — 1.24159 for Dy 
gives a reasonable minimum value to assume for the mag¬ 
netic moment difference, namely that with the nearest atomic 
Zeeman level. A good fit to the data is them obtained with 
ires = 16.17 for the 76G resonance. Similar quality fits can 
be obtained for other combinations of ires and dp such that 
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FIG. 12. (Upper panel pair) Scattering length and bound state en¬ 
ergies versus E^es for ires = 8.719. The dashed lines show the bare 
(uncoupled) energies and the solid red lines show the result of the 
coupled calculation. (Lower panel) Binding energy of the threshold 
bound state relative to its pole position at for three different val¬ 
ues of pole strength ires. The dashed line shows the ramping bare 
level that generates the coupled halo bound state. 




FIG. 13. (Upper Panel) Comparison of coupled channels model with 
the measured binding energies near the 76G resonance. Very simi¬ 
lar fits are obtained with different combinations of Jres ^ 1 and 5/i 
such that ires/(5/t) is constant. (Lower Panel) Similar fit of the cou¬ 
pled channels model to the measured binding energies near the 170G 
resonance. 


■Sres/(5/f) = /(a£) is constant. This determines abgA = 

2600 aoG and 2800 aoG for the two fits for the 76G res¬ 
onance. The calculated width of the resonance with Sres = 
16.17 is \Eres-E^,,,\/h = 47.11 MHz, or 27.1 G upon divid¬ 
ing hy 5fi = gUs- The fact that this value is consistent with 
the value of A = 31(6) G measured by examining the 3-body 
line shapes suggests that these values for ires and djJ. are plau¬ 
sible. A similar quality fit for the 170G resonance, shown in 
the figure, is obtained by taking the same ires values and as¬ 
suming a 5/r value 6 per cent larger. Thus, the two resonances 
have comparable pole strengths, not surprising since they have 
quite similar measured E versus B curves. 
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[1] T. Maier, H. Kadau, M. Schmitt, A. Griesmaier, and T. Pfau, 
Opt. Lett. 39,3138 (2014). 

[2] T. Maier, H. Kadau, M. Schmitt, M. Wenzel, 1. Ferrier-Barhut, 
T. Pfau, A. Frisch, S. Baier, K. Aikawa, L. Chomaz, M. J. Mark, 
F. Ferlaino, C. Makrides, E. Tiesinga, A. Petrov, and S. Ko- 
tochigova, (2015), 1506.05221. 

[3] U. Eismann, L. Khaykovich, S. Laurent, 1. Eerrier-Barhut, B. S. 
Rem, A. T. Grier, M. Delehaye, E. Chevy, C. Salomon, L.-C. 
Ha, and C. Chin, ArXiv.org (2015), 1505.04523. 

[4] A. Erisch, M. Mark, K. Aikawa, E. Ferlaino, J. L. Bohn, 
C. Makrides, A. Petrov, and S. Kotochigova, Nature 507 , 475 
(2014). 

[5] K. Baumann, N. Q. Burdick, M. Lu, and B. L. Lev, Phys. Rev. 
A 89, 020701 (2014). 















































15 


[6] M. Berninger, A. Zenesini, B. Huang, W. Harm, H.-C. Nagerl, 
F. Ferlaino, R. Grimm, R S. Julienne, and J. M. Hutson, Rhys. 
Rev. A 87, 032517 (2013). 

[7] K. Jachymski and R. S. Julienne, Rhys. Rev. A 88, 052701 
(2013). 

[8] S. Kotochigova and A. Retrov, Rhys. Chem. Chem. Rhys. 13, 
19165 (2011). 

[9] J. L. Bohn, M. Cavagnero, and C. Ticknor, New J. Rhys. 11 , 
055039 (2009). 

[10] B. Gao, Rhys. Rev. A 78, 012702 (2008). 

[11] C. Chin, R. Grimm, R. Julienne, and E. Tiesinga, Rev. Mod. 
Rhys. 82, 1225 (2010). 

[12] H.-W. Cho, D. J. McCarron, M. R. Koppinger, D. L. Jenkin, 


K. L. Butler, R. S. Julienne, C. L. Blackley, C. R. Le Sueur, J. M. 
Hutson, and S. L. Cornish, Rhys. Rev. A 87, 010703 (2013). 

[13] Y. Wang, R. S. Julienne, and C. H. Greene, ArXiv.org , 8094 
(2014), 1412.8094. 

[14] T. Nicholson, S. Blatt, B. Bloom, J. Williams, J. Thomsen, 
J. Ye, and R. S. Julienne, eprint arXiv: 1502.00026 (2015). 

[15] R. S. Julienne and J. M. Hutson, Rhys. Rev. A 89, 052715 
(2014). 

[16] F. H. Mies, E. Tiesinga, and R. S. Julienne, Rhys. Rev. A 61 , 
022721 (2000). 

[17] N. Nygaard, B. I. Schneider, and R. S. Julienne, Rhys. Rev. A 
73, 042705 (2006). 


